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I. SUMMARY 

The investigative efforts at Southwest Research Institute have centered 
on the development and execution of a comprehensive three-dimensional ray 
tracing code. This code contains a model of the Jovian magnetosphere which 
includes the 0-4 magnetic field description {Acuna and Ness, 1976) and an 
empirical plasma model. Decametric radiation (2-35MHz) has been the primary 
subject of the analyses but some studies in the whistler mode (< 10kHz) have 
been performed. The whistler mode analyses yielded the approximate dimensions 
of Jovian lightning source regions and the approximate energy dispersion. 

Prom the latter an estimate of the Jovian lightning emission rate was made. 
These results were previously reported in the Final Report submitted to the 
Jet Propulsion Laboratory in October, 1982. The decametric radiation (DAM) 
studies have produced numerous publishable results. 

The early work of the current contract in ray tracing of DAM did not 
incorporate a model Jovian ionosphere. Sources were assumed to be at the foot 
of the Io flux tube and only northern hemisphere sources were initially 
considered. Emission emanating in a hollow cone around the active field line 
was ray traced from sources separated by 10° in Jovian longitude. One of the 
most notable results was the observation that not only does the magnetic field 
model control the frequency extent of the emission but also the curvature and 
morphology of the arcs . The importance of displaying the Planetary Radio 
Astronomy data in a constant-lo longitude format were clearly demonstrated in 
this early work. The results were published in the Journal of Geophysical 
Research (Appendix 1). 

Isolated, high-curvature DAM arcs were modelled by varying the wave 
normal angle at the source until the ray tracing results matched the 
observations . Plots of wave normal angle versus frequency for different 
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Doppler shifts at the source were made which provide the first comprehensive 
three-dimensional description of the relationship of these variables. The 
results were compared with the theoretical work of Staelin (1981) and imply a 
maximum wave normal angle of about 80° for the high-curvature arcs. This work 
has been submitted to the Journal of Geophysical Research (Appendix 3) . 

Because the Jovian ionosphere was found to strongly influence the 
higher frequency radiation, an empirical expression was used to represent the 
ionosphere in the code. The ray tracing of Jovian DAM source points was 
repeated including the southern hemisphere sources. The results have, for the 
first time clearly delineated northern and southern hemisphere sources in the 
Voyager data. The southern hemisphere arcs are lower in frequency extent and 
have a narrower emission cone angle. These findings lend strong support to 
the theory that arc morphology is primarily a result of wave refraction in the 
magnetosphere rather than emission processes. Appendix 2 is preprint of a 
paper submitted to the Journal of Geophysical Research detailing these 
exciting results. 

The next step in reformatting the Voyager data led to a discovery of 
an lo-independent emission which is most probably the Jovian counterpart to 
terrestial auroral kilometric radiation (AKR) . By plotting the sub-io 
longitude versus the difference between the spacecraft and the sub-io 
longitude, an emission which has been named Jovian Auroral Decametric 
Radiation (JADR) has been identified. This emission appears to emanate only 
from high-latitude nightside auroral zone sources in a manner similiar to the 
terrestial AKR. These results have been submitted to the Journal of 
Geophysical Research and appear in preprint form in Appendix 4. 
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IX . INTRODUCTION 

This work represents culmination of a three year activity involving 
ray tracing studies of Jovian low frequency emisions. The program included 
the development and execution of a comprehensive three-dimensional ray tracing 
computer code for examination of model Jovian decametric (DAM) emission, a 
detailed description of this program has been submitted to the jet propulsion 
Laboratory as "Final Report" dated October, 1982. In the present report a 
brief description of the improvements to the computer code will be outlined 
and a presentation of the results of ray tracing of Jovian emissions will be 
presented in the form of summaries with details given in the appendices. The 
appendices contain copies of papers that have been accepted or submitted for 
publication to the journal of Geophysical Research. 
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III. TECHNICAL DISCUSSION 
A. previous Results 

As described previously in the "Pinal Report" dated October, 1982 
(hereafter referred to as FR1982), DAM emission from northern hemisphere 
source positions located at the feet of Io flux tubes separated by 10° in 
Jovian longitude have been ray traced. The results were compared to Voyager 1 
and 2 observations and have been published in the Journal of Geophysical 
Research . As a result of these efforts the importance of wave propagation 
effects in reproducing the observed Jovian DAM arcs became apparent. The 
magnetic field model and ionospheric effects were seen to be fundamentally 
important, in addition, the arcs were seen to be divided into high and low 
curvature classes . initial explanations for some of uhe differences between 
these classes were presented which included variations of initial wave normal 
'angle and local magnetic field strength at the source. Appendix I is a 
reprint of a paper published in the Journal of Geophysical Research entitled 
"Three Dimensional Ray Tracing of Jovian Radiation in the Low Frequency Range" 
which details the results of this original work. 

Initial studies of whistler mode emissions (l-3kHz) from lightning 
sources in the Jovian cloud tops were also made as discussed in FR1982. 
Whistlers were ray traced through the Jovian magnetosphere and into the Io 
torus in an effort to determine the extent of the source positions and the 
degree of energy dissipation from this emission. These results have been 
reported to the University of Iowa and are awaiting submission for 
publication. 
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B. Improvements to the Computer Code 

As detailed in FR1982, the ray tracing code is based on the Stix 
cold plasma formulation and solves the Haselgrove (1954) equations for the 
index of refraction. The method of solution is by Runge-Kutta integration 
with a modified Hamming predicator/corrector scheme. The ambient Jovian 
plasma density model is based on the Sentman-Goertz empirical results and the 
lo torus plasma density is a spline fit to the published contours of Warwick 
et al. (1979). The magnetic field model used is the 0-4 model of Acuna and 
Ness (1976). Now included in the code is an ionospheric model based on the 
empirical nonlinear fitting of Hashimoto and Goldstein, (1983). 

Atreya and Donahue, 1976 have derived a model Jovian ionosphere 
based on observations. Hashimoto and Goldstein, (1983) have used a nonlinear 
least squares fit to the model of Atreya and Donahue producing the following 
analytical result! 

N x “1 .15x105 exp((1417-h/386.3)/(1+(991/(h-276)) 4 * 31 ), 
where N x (cm~ 3 ) is the ionospheric electron density, h is the altitude in km 
above the cloud tops and N x =0 for h < 276 km. The magnetospheric model and 
the ionospheric models have been matched at the top of the ionosphere so that 
a smooth transition occurs. With the addition of the ionosphere it was 
necessary to ray trace DAM at all frequencies whose source positions reached 
the ionosphere. 

In order to display the results of the ray tracing code a plot 
program was developed on an HP1000 computer which provided high quality ink 
bed plots. The wave frequency versus Jovian system III longitude was 
displayed on these plots with sources for each sub-Io longitude clearly 
indicated by color symbols . This plotting technique allowed easy comparison 
of the model arcs with Voyager 1 and 2 observations. 
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C. New Results 

As a result of the identification of the Jovian ionosphere as an 
important source of wave refraction at certain DAM frequencies, it was 
necessary to ray ii‘ace all the northern hemisphere DAM source positions again. 
The results of Staelin, (1981) suggested a constant initial wave normal angle 
of 85° at the source (instead of the original 90° used in the previous study 
[cf. Appendix 13). Figure 1 (A-D) shows the resulting model emission on a 
frequency-versus-system Ill-longitude scale for arcs within the longitude 
ranges of 10°-90°, 100°-180°, 190°-270°, and 280°-360° respectively. Figure 2 
is in the same format for the source points of the southern hemisphere. It is 
immediately apparent that the model southern hemisphere source emissions have 
a lower frequency cutoff than the northern hemisphere sources, and emission 
cone edges are closer together in longitude. In order to directly compare the 
model ray tracing results with actual observations it was decided to re-format 
the Voyager 1 and 2 Planetary Radio Astronomy (PRA) data by sorting it with a 
fixed sub-Io longitude . This allowed direct comparison with the ray tracing 
results. The reformatted data was obtained from Marshall Space Flight Center 
as a result of the efforts of J. L, Green and N. F. Six, supported by the 
Jupiter Data Analysis program. Comparison of the model and observations was 
very good for some cases and allowed the first delineation of Jovian DAM 
emissions into southern and northern hemisphere sources. A detailed 
discussion of the comparison is included in this report as Appendix 2, which 
is a preprint of a paper entitled "Identification of Decametric Radiation From 
Sources in the southern Hemisphere of Jupiter" submitted for publication to 
the Journal of Geophysical Research. 
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Careful study of the dynamic frequency- time voyager spectrograms reveals 
a multitude of arc morphologies with arcs varying in frequency extent and in 
curvature. Included is a class of arc with large frequency extent and 
curvature which appears often as intense, and somewhat isolated from the 
surrounding arcs. By assuming a wave normal angle which varied with frequency 
at the source point we were able to match the curvature of these observed 
high-curvature arcs with model ray tracing results. The wave normal angles 
chosen were consistent with relativistically Doppler-shifted electron emission 
in RX mode {staelin, 1981). The electron energy was found to be about 6 kev 
which is in agreement with electron energies in the vicinity of the Io torus 
which are hypothesized to be accelerated and produce the DAM. This work 
represents the first comprehensive three-dimensional ray tracing analysis of 
varying the wave normal angle, previous studies considered only a 
two-dimensional analysis and consequently could not adequately reproduce 
observed arcs. The results of the present analysis are detailed in Appendix 3 
which is a preprint of a paper entitled "Jovian Decametric Arcs; An Estimate 
of the Required Wave Normal Angles From Three-Dimensional Ray Tracing" , 
submitted to the Journal of Geophysical Research . 

The reformatting of the Voyager Pra data has allowed rapid analysis of 
both io-dependent and lo-independent emissions . Carrying the process one step 
further, s. Gulkis plotted the data using the sub-lo longitude as the 
ordinate and the difference between the spacecraft longitude and the sub-lo 
longitude as the coordinate. This technique immediately allowed the 
identification of an lo-independent emission which has been named Jovian 
Auroral Decametric Radiation (JADR) . This radiation is seen only on the 
nightside and appears to have a source in the high-latitude auroral zone. 

These facts indicate the emission may be the Jovian counterpart to the 
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intensely studied terrestial AKR emissions. Appendix 4 i3 a preprint of a 
j&ipar entitled "identification of a Night-Side Component of DAM as a Jovian 
Counterpart to AKR" submitted to the Journal of Geophysical Research which 
details the exciting discovery. 
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IV, CONCLUSIONS 

Results of three-dimensional ray tracing of Jovian low frequency 
radiation have provided important information regarding l&ija propagation of DAM 
emissions. Direct comparison of model ray tracing results with observed arcs 
has indicated that many arcs can be adequately explained by propagation from 
hollow emission cones with sources at the feet of Io flux tubes. Emission 
occurs at frequencies Blightly Doppler-shifted above the RX cutoff frequency, 
with wave normal angles close to 90°. In addition, the effect of the 
magnetic field in determining the curvature and frequency extent of the arcs 
is fundamental (cf. Appendix 1). 

After the introduction of a Jovian ionospheric model new ray tracing 
results have yielded improved model arcs for the northern hemisphere and the 
first three-dimensional model arcs with source positions in the southern 
hemisphere (Appendix 2) . The ray tracing has allowed the first delineation of 
northern and southern hemisphere emissions of Voyager data. The importance of 
re-formatting the Voyager data 3uch that the sub-Io longitude can be held 
fixed has been found to be of fundamental importance in determining 
Io-dependent spectral features. 

Certain high-curvature arc signatures can be modelled by allowing the 
wave normal angle to vary yielding a ¥(f) relationship that is unique to the 
high-curvature arcs. One explanation for the smaller wave normal angles 
necessary to model certain high-curvature arcs is relativistically 
Doppler-shifted electron gyro-emission (Appendix 3). For electrons with 
energies of about 10 kev (reasonable for the source particles) the wave normal 
angles resulting from theoretical Doppler shifted emission are quite close to 
the values obtained from ray tracing needed to match the observed arcs . The 
implication is that certain high-curvature arcs result from the gyro-emission 
of beamed electrons with energies of a few kev. 
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Reformatting the Voyager PRA data has also allowed the identification 
of an Io-independent emission which is the probable couterpart to the 
terrestial AKR. Jovian Auroral Decametric Radiation may provide new insights 
to the emission mechanism responsible for the production of AKR, as well as 
provide valuable information about the structure of the Jovian magnetosphere. 


IX 
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FIGURE CAPTIONS 

FIGURE Is Frequency versus system ill '65 longitude showing the model arcs 
at 10° intervals of constant Io CML. The Io CML is labelled on 
each curve. The curves represent the intersection of one "side" 
of the emission cone with a radial shell at 150 Rj . The source 
points for the arcs shown lie on io flux tubes in the northern 
hemisphere with longitudes in the range: A) 10° to 90°; 

B) 100° to 180 ° t C) 190° to 270°; and D) 280° to 360°, 

FIGURE 2: Same as Figure 1 except the source points for the arcs shown lie 

on Io flux tubes in the SOUTHERN hemisphere with longitudes in 
the range: A) 10° to 90°; B) 100° to 180°? C) 190° to 270°; and 

D) 280° to 360°. 
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A. Appendix 1i "Three-Dimensional Ray Tracing of 
Jovian Radiation in the Low Frequency Range" 
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Three-Dimensional Ray Tracing of the Jovian 
Magnetosphere in the Low-Frequency Range 

J. D. MENtETTi,* J. L. Green. 1 2 S. Gijlkis, 3 and F, Six 4 

Three-dimensional ray tracing of the Jovian DAM emission has been performed utilizing the 0*4 
magnetic held model (Acuna and Ness. 1979) and a realistic plasma model. Minimal assumptions about 
the emission mechanism have been made that include radiation in the right-hand extraordinary mode, 
propagating nearly perpendicular to the Held line at source points located just above the RX cutoff 
frequency along lo flux lubes. We have performed ray tracing In the frequency range from 2-35 MHz 
from successive lo flux tubes separated by ten degrees of central meridian longitude for u full 
circumference of nonhern hemisphere sources. Our results show unusual complexity In model arc 
spectra that Is displayed in a constant lo phase format with many similarities to the Voyager PRA data, 
The results suggest much of the variation In observed DAM spectral features is a result of propagation 
effects rather than emission process differences. We feel our work represents a comprehensive ray 
(racing Investigation of DAM arc emission and presentation of the results In a unique format. 


Introduction 

One of the most characteristic features of the Jovian 
decametric radio emissions is the presence of arc*like struc- 
tures in frequency-time spectrograms. The arc-like struc- 
tures were first observed by the Voyager planetary radio 
astronomy experiment (PRA) [Warwick et at,. 1979] and 
later found in ground based observations [Bolschot et at., 
1981). The arcs consist of narrow band emissions that drift in 
frequency in a repeatable and highly organized manner as 
shown in Plate 1. The frequency range of the arc emission is 
from 40 MHz down to several MHz or lower, but all arcs do 
not span this full frequency range. Several different kinds of 
arcs have been recognized in the data, depending on their 
curvatures, which either open toward increasing time (ver- 
tex early) or open toward decreasing time (vertex late). The 
arcs often appear as a closely spaced series referred to as 
"nested arcs." Detailed descriptions of the arcs can be 
found elsewhere fi.o., Warwick et at., !979; Leblanc, 1981; 
Bolschot and Aubtcr, 1981, Carr et al„ 1983). 

Various source mechanisms have been proposed to ex- 
plain the decametric arcs. A general feature of most theories 
is that the emission occurs along conical sheets that are 
swept past the observer as the planet rotates && first suggest- 
ed by Dulk (1965). Many investigators attribute the emission 
to precipitating electrons along the lo magnetic flux tube that 
radiate near the gyrofrequency [ Coldrelch and Lyndcn-Bcll. 
1969; Smith, 1976; Goldstein and Gocrtz, 1983). Goldstein 
and Theiman [1981J, Pearce (1581), Staelln f 1981 ). Neu- 
baner [1980] have shown that arc-like structures arc pro- 
duced by the conical sheet model although no simple choice 
of parameters has been found which reproduces the data in 
detail. In particular, the cone emission model does not 
explain the "nested arcs." These have been attributed to 
multiple reflections of Alf''?.n waves [Garnett and Gocrtz, 
1981), interference effects [Bolschot and Aubicr, 1981 ; Leva- 
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cheux et al„ 1981). and local magnetic anomalies [Warwick, 
1981). 

Hashimoto and Goldstein f 1982] used a three-dimensional 
ray tracing program to investigate the observed occurrence 
probability maps of the lo-dependent Jovian decametric 
radiation, They also Investigated conditions for wave gener- 
ation in the Jovian ionosphere. Goldstein and Theiman 
[1981] performed an analytical study of refraction in the 
Jovian magnetosphere and ionosphere. They demonstrated 
that an extraordinary wave. Initially excited nearly orthogo- 
nal to the magnetic field, is significantly influenced by the 
local plasma density and that arc-like features (hat resemble 
those observed are produced by the geometry of the model, 
lit this paper, we extend the work of Hashimoto and 
Goldstein [1983) and Goldstein and Theiman (1981) by 
studying the arc structures using a three-dimensional ray 
tracing program. 

The purpose of this paper is to perform three-dimensional 
ray tracing of DAM radiation from source positions along an 
lo flux tube at frequencies greater than but nearly equal to 
the right-hand extraordinary (RX) mode cutoff frequency 
/rx a IA + «/*/ 2) : + where f„ is gyrofrequency and 
Jp is plasma frequency. We have incorporated a realistic 
Jovian plasma mode! [Sentman and Goertz, 1978], and the 
0-4 magnetic field model [Acuna and Ness, 1976], Jovian 
DAM arcs are produced with minimal assumptions about the 
source mechanism. 

Rav Tracing Equations and 
Magnetospheric Model 

Ray Tracing 

To calculate propagation ray paths in the Jovian magneto- 
sphere. we have used a three-dimensional ray tracing pro- 
gram. The computer code is based on the Stix [1962] coid 
plasma formulation of the index of refraction and Hasel- 
grove's [1955] set of first order differential equations amena- 
ble to numerical solution. The program was initially devel- 
oped by Shawhan [1966] for two-dimensional ray tracing in 
the terrestial ionosphere and magnetosphere. For the pur- 
pose of this three-dimensional ray tracing study, we have 
incorporated additional equations that allow for variations in 
the azimuthal direction. This program has been extensively 
tested and reduces to the two dimensional results of Green et 
al. [1977], Green and Garnett [1980], and Lecacheaax 
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[1981,. Chosen as inputs are a source position on radial 
distance, longitude, and latitude) from which waves at a 
specific frequency / are to be launched at a specific wave 
normal angle Ji. with respect to the magnetic held Different 
frequencies originating from source points on ihe same held 
line or on adjacent held lines would generate a set of nested 
emission cones Depending on (he geemetry of these cones, 
imposed by the plasma and by the B held, the spacecraft 
would detect frequencies that decrease or increase with time 
as these cones sweep over the detector 

At the source region, the program calculates an index of 
refraction surface tor RX mode waves, based on held and 
on background plasma parameters Next, the program takes 
an increment*' step in the direction perpendicular to 'he 
index of refraction surface. : e . in the direction of the group 
velocity of energy .'„w. and then determines the coordinates 
of this new point on the ray path Then, another index of 
refraction surface is calculated through this new point and 
the steps are repeated The direction of each ray changes 
according to Snell's law as it travels through the magnetized 
plasma To generate an emission cone, thirty-six rays are 
launched, one every 10” around B at the specified wave 
normal angle and frequency. 

B field and Plasma Models 

The ray paths calculated by the procedure described 
above are dependent on the magnetosphenc model chosen. 
The description of the Jovian magnetic fHd utilized in these 
calculations (the 0-4 model) is that published by Acuna and 
Ness [ I976|; the description of the background plasma is that 
published by Seniman and Goertz (1978), and the descrip- 
tion the lo torus is that published by Warwick et al. (I979|. 
Figure I is an abbreviation of this combined "Jovian magne- 
tosphere" showing contours of plasma frequency in the 200“ 
-2C” mendian plane Suspected source regions in the north- 


ern and southern hemispheres are shaded Sphencal har- 
monic expansion coefficients are used to calculate the mag- 
nitude of the field B n . the c smponents Br. B and B t and 
their derivatives The plasma density along the ray path, and 
the spatial gradients are calculated from -mpirtcal fits 

As Hashimolo and Goldstein |I983| have pointed out. 
large plasma densities can produce important refractory 
effects at the higher frequencies when source points near the 
foot of the lo flux tube extend down into the ionosphere We 
have incorporated the ionospheric model introduced by 
Hashimolo and Goldstein to determine the upper cutoff 
frequency associated with each lo flux tube investigated in 
this study, but we have not incorporated an ionospheric 
model in our ray tracing studies We have avoided ray 
tracing at frequencies where ionospheric influence is impor- 
tant We became aware of the importance of the ionosphere 
in affecting propagation at certain frequencies only after our 
extensive ray tracing investigation had aJmost been complci 
ed. Tests at several individual lo flux tubes have indicated 
that incorporating an ionospheric model that continuously 
joins with a magnetosphenc model will modify our results 
only slightly at frequencies whose source points lie outside 
the influence of the ionosphere 

Computations 

In addition to the magnetosphenc plasma model discussed 
in the previous section, some assumptions regarding source 
location, initial launch angle, and mode of propagation are 
required as inputs to the ray tracing studies. Rather than 
carrying out a parameter study of these vanables. we have 
chosen a set of initial conditions and studied the propagation 
over a full range of observer-lo-Jupiter longitude phase 
space. The results of our studies are intended to show the 
trends in the radiation spectra introduced by the complex 
geometry of the problem Indeed, as we will discuss below. 
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SYSTEM III ( 1969) LONGITUDE 

Fi| J. Calculated radiation sheen for two northern hemisphere source* i2 MHz and 10 MHr) and lor two southern 
hemisphere source* (2 MHz and 20 MHzi The radial shell al 100 Jupiier radii 1 100 R,) is shown in mercaior projection 
All ray* were launched al 90* with respect lo .he B held 


our studies show that exceedingly complex spectra can be 
produced by propagation effects alone, without introducing 
complexity into the emission mechanism itself. 

The initial condition* we used are a* follow*: 

I The source region is at the fool of the magnetic L shell 
that intercept* lo 

2. The frequency of emission in the source region is 
close to the gyrofrequency and greater than the RX cutoff 
frequency. 

3. The emitted radiation is in the right-hand polarized 
extraordinary mode 

4 The radiation is emitted in a hollow cone pattern about 
the magnetic held direction, with the emission angle close to 
90° (with respect to the magnetic held). 

No specihc emission mechanisms are implied by our as- 
sumption: rather, they are chosen to be consistent with (he 
available data (i.e., Carr el al., 19831 and to at least some of 
the theories of radio emissions (i.e . Dulk. 1965. Seubauer. 
1980. Goenz. 1980 , Gurnelt anj Goeriz. 1981, Goldurtn and 
Goem. 1983). 

An entire senes of ray trajectones are computed for lo 
positioned every 10* of Jupiter longitude Typically, lo's 
position is chosen and a magnetic line of force is calculated 
from lo to the foot of the flux tube near Jupiter's cloud tops. 
Source points are located at a set of frequencies, typically 2. 
5. 10. IS. 20, 25 . 30. and 35 MHz from the condition ///, - 
1 .005 al the origin. This latter ratio is chosen to be consistent 
with slightly Doppler-shifted emission al large wave normal 
angles. Gyroemission is not assumed. The RX cutoff fre- 
quency is calculated and compared with /. Under the 
condition / > / HX . wave propagation is possible and a senes 
of rays is traced at 36 evenly spaced azimuthal directions 
around the source point. We have also accounted for the 
oblateness of Jupiter in our determination of /*x The radial 
distance to the cloud tops, fl(km) is give" by fi * 71300(1 - 
cos : (colatitude»/l5.4). Figure I shows on a mercator oroiec- 
non. the intersection of several families of ray trajectories 
with a sphere of radius 100 R, that surrounds Jupiter Shown 
in the figure are a 2- and 10-MHz family that onginate in the 
northern hemisphere, and a 2- and 20- MHz family that 
onginate in the southern hemisphere. 


In all. more than 1(8) hours of computer time (CPU) on a 
CYBER 172 was required for the results reported in this 
study Typically. 23 mm frequency for each lo position was 
required (o track the ray from the source point out to 150 R,. 

Results 

The principal results of our computations are shown in 
Figures 3 a and 3 b The figure shows the frequencies that 
would intersect a spacecraft were it located at R ■ 150 R/ 
and at a latitude of 3 2* The curves are parameterized by the 
central mendian longitude facing lo at the time of emission. 
At this distance from the source location, the parallax 
between the center of the planet and (he source is small, and 
the results apply generally to any observer on this latitude 
and at a greater distance from Jupiter. The two figures 
correspond to the two azimuth positions that intersect the 
observer from every source point (i.e.. in general, emission 
from a cone may intersect the observer in two azimuth 
directions) The solid curves in Figures 3u and 36 are not in 
general what a fixed observer would see because lo is not 
stationary in the Jovian-fixed longitude system. Hence a 
stationary observer sees lo moving from longitude to longi- 
tude. 

Discussion 

In Figures 3o and 3 b are revealed the wide range of 
structure in the curves resulting from varying the lo flux tube 
longitude of the source points. Indicated by small horizontal 
brackets on Figures 3a and 3 b are the upper cutoff frequen- 
cies for RX mode emission at source point frequencies 
slightly above the RX cutoff. These cutoffs indicate the 
maximum frequency that is expected to be emitted along a 
particular lo flux tube if that emission is slightly above fn%. 
No DAM emission will escape the ionosphere if the source is 
at an altitude less than the altitude of the density maximum 
in the ionosphere, We have so far not performed ray tracing 
at frequencies that are determined to have sources deep in 
the ionosphere but have used an ionospheric model only to 
determine what frequencies will be influenced and the cutoff 
frequency. Oftentimes, however, one does not observe 
emission at higher frequencies because those frequencies are 
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Fig \u Frequency versus system III 65 longitude shotting the model arcs at 10 intervals of constant lo CML 
(constant lo format i The loCMl. is labeled on each curve These curses represent the intersection of one "side ' of the 
emission cone with a radial shell at 150 H , The horizontal brackets indicate the maximum frequency that can be emitted 
under the conditions specified in the text Also indicated are the location of the lo-B source and a modified sketch of 
part of Figure 5 |S [Cun n ul , IVKi | indicating the typical signature seen on a dynamic tnonconstant lo) format 


refracted (without ionospheric influence) and do not inter- 
cept the spacecraft The upper cutoff frequency indicated on 
Figures 3« and 36 are thus a result of the complex interaction 
of the magnetic field structure, ionosphere, magnetosphere, 
and oblateness of Jupiter 

In Figures 3 a and 3 b for lo central meridian longitude 
(CML) in the range from approximately 80‘ to 240 the 
curves are rather straight, indicating thai multiple frequen- 
cies may be observed al a given spacecraft CML By "lo 
CML" we mean the sub-Io longitude For lo CML in the 
range from 360' through 70 c the curves have greater slope, 
and some show more curvature with a lower but increasing 
upper cutoff frequency In the range of lo CML from 250‘ 
through 350‘. the upper cutoff frequencies are lower with 
complex structure shown in Figure 3 b when the spacecraft 
CML is greater than the CML of lo, It should be noted that 
while lo may be located at a given CML. the foot of the lo 
flux tube is not. in general, at the same CML. We believe 
much of the changing morphology seen in the curves show n 
in Figures 3a and 3 b is due to the azimuthal asymmetry of 
the 0-4 magnetic field model, which changes the source 
point field strength and latitude in a circumference of Jupi- 
ter. 

Also shown superimposed on Figures 3o and 3 b are 
modified sketches of part of Figure 7.15 from Carr el al. 
(1983). These sketches indicate the location and typical 
signatures of the lo-dependent dynamic (non-lo-fixedl spec- 
tra. The lo-B source, for instance, generally displays fairly 
low curvature arcs with smoothly falling upper cutoff fre- 


quencies in general agreement with the model arcs of Figure 
3 a The lo-A source typically displays arcs of more curva- 
ture than lo-B source and also displays a lower cutoff 
frequency with increasing spacecraft CML The model arcs 
of Figure 3 b in the lo-A source region would map to lower 
curvature arcs on a dynamic tnon-lo-fixed) plot, but the 
upper cutoff frequency does fall off as observed For the lo-C 
source the overlap between the observed lo phase and the 
model lo phase is not as good as for the lo-A and B sources 
(cf Table 7 4 of Carr et al. |I983|). but the Io-C source 
region includes both left-hand and right-hand polarization 
emissions, which indicates possibly some southern hemi- 
sphere sources that were not considered in the production of 
Figure 3 The lo-C source region dynamic spectra show large 
curvature with rapidly dropping upper cutoff frequency in a 
region of Figure 3 h where the model curves show the 
greatest curvature and a sharp drop in the cutoff frequency . 

While the exact arc shape as noted on Figures 3 a and 3 b 
drawn in a constant-lo format should not be directly com- 
pared with the dynamic spectra that are drawn in a noncon- 
stant lo format, it is illustrative to note the changes in lo- 
dependent arc morphology as one moves to the different lo- 
A. -B. and -C source regions on the constant-lo format on 
Figures 3 a and 3 b It is tempting to infer that the morphologi- 
cal changes in lo-dependent arc structure as seen on the 
dynamic spectra can be largely explained by propagation 
effects rather than by varying emission processes 

Plotting the Voyager PRA data in a constant lo format can 
be a valuable tool to a better understanding of the emission 
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Fig. 3 b Some as Figure 3a for the other "side'' of the emission cone The lo-A and lo-C source locations and typical 
dynamic arc signatures are also indicated [Carr tt at . 19831 


processes of the lo-depcndent sources [cf. Green, I983|. We 
have avoided replotting Figure 3 in a dynamic -spectra format 
for direct comparison with the observed Voyager dynamic 
spectra (although this has been done for specific examples 
with reasonably good agreement), because we feel the value 
of Figure 3 is in displaying the major features of lo-depen- 
dent emission. The correlation found between Figures 3a 
and 3b and the lo-A, B. and -C source regions supports the 
initial assumptions pertaining to the emission process, which 
include RX mode emission almost perpendicular to the 
magnetic field from source points just above the KX cutoff 
We are currently in the process of comparing refoi-matted 
Voyager PRA data with Figure 3 in an effort to understand 
better the emission process(es) and to explain on a finer 
scale any differences noted We recognize these differences 
may be due to varying wave normal angle | Goldstein anti 
Thieman, 1981 1 and/or discrepancies in the magnetic or 
plasma field models. In addition, the ionosphere is row 


known to be of importance at the higher frequencies at some 
CML values. 

Conclusions 

We conclude from our results that with minimal assump- 
tions pertaining to an emission mechanism our ray tracing 
results produce a complex spectrum of Io-dependent DAM 
arcs as shown in Figures 3a and 3b These few assumptions 
together with the incorporation of the 0-4 magnetic field 
model [Acuna and Mess, 1979) and a plasma model due to 
Sentman and Goertz [ 1 9781 have resulted in model DAM 
arcs which display different curvatures and upper cutoff 
frequencies. 

When the source points are located on Io flux tubes that 
are at Jovian system III longitudes ranging from approxi- 
mately 90*-230*. where the 0-4 model yields high surface 
magnetic field strengths. Figures }a and 3b indicate that the 
resultant DAM arcs have low curvature and upper cutoff 
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frequencies extending beyond 35 MHz. On the other hand, 
when the fool of the lo flux tubes lie in the range of system 
III longitudes extending from 240'-80 # . where the 0-4 sur- 
face field strengths reach a minimum, the resultant arcs have 
lower upper cutoff frequency and a higher curvature. It is 
apparent that the B field structure that produces source 
points at the foot of the Io flux tube ranging in latitude from 
4S°-70° and in altitude from 0 to >2 Rj around the circumfer- 
ence of Jupiter produces substantial differences in the result- 
ant DAM arcs. 

We have also performed preliminary ray tracing studies, 
including the ionospheric model introduced by Hashiitioto 
and Goldstein [1982] and find that the ionosphere can have 
great influence on emission ray path at the higher frequen- 
cies of emission along some lo flux tubes. 

Our studies indicate many areas of future attention. We 
recognize the importance of including the ionosphere in our 
ray tracing and have begun reproducing our results including 
an ionosphere for both northern and southern hemisphere 
sources. In addition, we are investigating the effects of 
varying the wave normal angle and Doppler shift at the 
source point. We hope to explain better the higher curvature 
arcs often observed in the Voyager PRA data. 

Because of our apparent success at reproducing, at least 
qualitatively, some of the Voyager lo-dependcnt DAM ob- 
servations. we recognize the importance of reformatting the 
Voyager PRA data in a manner similar lo Figures 3 a and 
3 b — holding the lo longitude constant [cf. Green, 1983J. 
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frequencies extending beyond 35 MHz On Ihe other hand 
when the foot of the lo flux tubes lie in the range of system 
III longitudes extending from 240“-80' where the 0-4 sur- 
face held strengths reach a minimum, the resultant arcs have 
lower upper cutoff frequency and a higher curvature It is 
apparent that Ihe B held structure that produces source 
points at the foot of the lo flux tube ranging in latitude from 
4i*-70‘ and in altitude from 0 to >2 Rj around the circumfer- 
ence of Jupiter produces substantial differences in the result- 
ant DAM arcs 

We have also performed preliminary ray tracing studies, 
including the ionospheric model introduced by Hashimrsto 
and Goldstein 1 1*4821 and find that the ionosphere can have 
great influence on emission ray path at the higher frequen- 
cies of emission along some lo flux tubes 
Our studies indicate many areas of future attention We 
recognize the importance of including the ionosphere in our 
ray tracing and have begun reproducing our results including 
an ionosphere for both northern and southern hemisphere 
sources In addition, we are investigating the effects of 
varying the wave normal angle and Doppler shift at the 
source point We hope to explain better the higher curvature 
arcs often observed in the Voyager PR A data 
Because of our apparent success at reproducing, at least 
qualitatively, some of the Voyager lo-dependent DAM ob- 
servations. we recognize the importance of reformatting the 
Voyager PRA data in a manner similar to Figures 3 u and 
3b — holding the lo longitude constant |cf Green. 19831 
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ABSTRACT 


A technique has been developed which aids the identification of 
io-dependent decametric radiation originating from the southern hemisphere of 
Jupiter. This technique compares the results of model ray tracing 
calculations with the Planetary Radio Astronomy (PRA) observations. All 
Voyager 1 or 2 PRA observations are sorted into bins (±3° wide) centered on 
a specific lo cental meridian longitude. When the data is plotted (as a 
frequency-longitude spectrogram) in this coordinate system, io-dependent 
features can be identified and used for direct comparison with model ray 
tracing calculations. The ray tracing calculations are done in a model Jovian 
magnetosphere where it is assumed that the decametric emissions are generated 
in the RX mode from low altitude source regions along the instantaneous lo 
flux tube, in this study, we compare the observations taken for four specific 
lo longitudes (150°, 160°, 170°, and 260°) with the corresponding model ray 
tracing. The Voyager observations and the model ray tracing calculations, 
when combined, clearly distinguish the northern and southern hemisphere 
emission features. This method allows, for the first time, Io-dependent DAM 
emissions generated in the southern and northern hemisphere to be delineated. 



INTRODUCTION 


The complete spectrum of Jovian decametric (DAM) radiation was first 
observed by the Voyager Planetary Radio Astronomy (PRA) experiment (Warwick et 
al., 1979). One of the most striking features of the Jovian DAM amissions 
seen in frequency-time spectrograms, are "arc-like" bands, The frequency 
range of the arc emissions is from 1 to almost 40 MHz, but not all arcs span 
this complete frequency range. Several different types of arcs have been 
recognized in the data, depending on their curvatures which either open toward 
increasing time (vertex early) or open toward decreasing time (vertex late). 

In general, an arc may appear either as an isolated feature or there may be a 
series of closely spaced arcs (called "nested arcs"). Detailed descriptions 
of the Jovian decametric arcs can be found elsewhere (ie,, Warwick et al., 
1979; Leblanc, 1981; Boisehot and Aubier, 1981; Carr, Desch, and Alexander, 
1983) . 

A general feature of most propagation theories of the DAM arc producing 
mechanism is that the emission occurs along conical sheets that are swept past 
the observer as the planet rotates as first suggested by Dulk (1967). Many 
investigators attribute the emission to precipitating electrons along the io 
magnetic flux tube which radiate near the gyrof requency (Goldreich and 
Lynden-Bell (1969); Smith (1976); Goldstein and Goertz, 1983). Goldstein and 
Thieman (1981); Pearce (1981); Staelin (1981); Neubauer, (1980) have shown 
that arc-like structures are produced by a conical sheet model although no 
simple choice of parameters has been found which reproduces the data in 
detail, in particular, these simple emission cone models do not explain the 
"nested arcs". However, nested arcs have been attributed to multiple 
reflections of Alfven waves (Gurnett and Goertz, 1981), interference effects 
(Boisehot and Aubier, 1981; Lecacheux, 1981), and local magnetic anomalies 



(Warwick, 1981), As shown by Goldstein and Thieman (1981) and Menietti et al. 
(1984), the arc curvature which varios on a frequency-time spectrogram from 
almost zero to quite dramatic loops (great arcs) may result from variations of 
the wave normal angle as a function of frequency in the source region. 
Hashimoto and Goldstein (1983) have suggested that the lo pha.'}* asymmetry of 
the DAM arcs may be due to the difference in propagation time of Alfven waves 
reflected from both the northern and southern hemisphere Jovian ionosphere. 

The purpose of this paper is to compare three-dimensional ray tracing of 
Jovian DAM emissions, from the southern and northern hemisphere, with the PRA 
measurements in order to determine the frequency-time signatures of southern 
hemisphere arc3. This paper complements the northern hemisphere ray tracing 
of Menietti et. al. (1984). The observations will be transformed into an 
Xo-jupiter stationary coosdinate system (see Green, 1984). This coordinate 
system allows for the lo-dependent emission features observed by the PRA 
experiment to be directly compared with the model ray tracing results. We 
have chosen source regions at four sub-Io longitudes. The agreement between 
the Voyager observations and the model ray tracing is good and allows the 
separation of the northern and southern hemisphere emission features. With 
this technique the southern and northern hemisphere decametric sources have, 
we believe, for the first time been delineated. The good agreement between 
the observations and the ray tracing results provides strong evidence that the 
basic assumptions of the model are correct. Refinements are necessary, 
however, to reproduce exact details of the emissions. 



THE RAY TRACING PROGRAM 


The ray tracing program used for this study is discussed in some detail 
in Menietti et al. <1984). Briefly, the Jovian magnetosphere used is the 0-4 
magnetic field model of Acuna and Ness (1976), and the Jovian background 
magnetospheric plasma is given by the model proposed by Sentman. and Goertz 
(1979). The io torus plasma density is a spline fit to the published contours 
of Warwick et al, (1979). in addition, the computer code contains an 
empirical model of the Jovian ionosphere based on the results of Hashimoto and 
Goldstein (1983). The program is based on a closed set of first order 
differential equations in spherical polar coordinates specialized by Shawhan 
(1966) (in two dimensions) and by Menietti et al. (1984) (in three dimensions) 
that describe the path of rays in a magnetized plasma. The expressions for 
the phase index of refraction and its derivatives are in the cold plasma 
formulation of Stix (1962). 

Qualitatively, the program calculates an index of refraction surface and 
group velocity (magnitude and direction) for a predetermined plasma wave mode 
in the model magnetosphere at some starting point. The program then takes an 
incremental step in the direction of the group velocity or energy flow and 
determines the coordinates of a new point on the raypath. Another index of 
refraction surface is calculated at this new point and the steps are repeated. 
The program used here has been extensively tested and reduces to the 
two-dimensional results of Green and Gurnett (1980) and Lecacheux (1981). 

MODEL ASSUMPTIONS 

The following assumptions regarding source location, initial launch angle 
and mode of propagation are required as inputs to the ray tracing program: 



1 . The source region i3 at the foot of the magnetic L-shell which 
intercepts lo; 

2. The frequency of emission in the source region is close to the 
gyrofrequency and greater than the RX cutoff frequency; 

3. The emitted radiation is in the right-hand polarized extraordinary 
mode; and 

4. The radiation is emitted in a hollow cone pattern about the magnetic 
field direction. The emission angle with respect to the magnetic 
field was set equal to 85° for northern hemisphere sources and to 
both 85° and 95° for southern hemisphere sources. 

No specific emission mechanisms are implied by our assumptions; rather, 
they are chosen to be cbnsistent with the available data (i.e. Carr et al., 
1983) and at least some of the theories of radio emissions (i.e. Dulk, 1967; 
Neubauer, 1980, Goertz, 1980; Gurnett and Goertz, 1981; Goldstein and Goertz, 
1983) . 

We have computed ray trajectories for sources located at sub-lo 
longitudes located every 10° of Jovian longitude, for both northern and 
southern hemisphere sources. Source points are located along an Io flux tube 
at a set of frequencies, typically 2, 5, 10, 15, 20, 25, 30, and 35 MHz from 
the condition f/f g =»1 *01 at the the origin. This latter ratio is chosen to be 
consistent with slightly Doppler-shifted emission at a wave normal angle, Y 0 i 
of 85° (cf. Staelin, 1981), Gyroemission is not assumed. Under the condition 
that the wave frequency is greater than the RX cutoff frequency, wave 
propagation is possible and a series of rays, for a constant frequency, is 
traced at 36 evenly spaced azimuthal directions around the magnetic field at 
the chosen source location. In the determination of the location of the RX 



cutoff frequency, the oblateness of Jupiter is accounted for; the radial 
distance to the cloud tops is given by R ■ 71300 (1~cos 2 {colatitude)/15.4) . 

The intersection of the resulting generated emission cones, for all 
frequencies, at 150 Rj and at a latitude of 5° determines the model decametric 
emission spectrum. It is this model decametric spectrum which is compared 
with the observations. For the purposes of this study we have chosen only 
four sub-Io longitudes to compare directly with observations which we feel 
clearly depict those emissions from southern hemisphere sources. 

COMPARISON OF OBSERVATIONS AND RAY TRACING RESULTS 

Only PRA observations made from February 21 to March 18, 1979 by Voyager 
1 and from June 26 to July 23, 1979 by Voyager 2 are used in this study. 

These observations are approximately centered around the encounter period and 
comprise enough data for a complete frequency-longitude study of lo-dependent 
DAM emissions to be made, observations at radial distances of less than 30 Rj 
wore excluded in order to avoid the changing emission pattern due to 
inner-magnetospheric propagation effects. The frequency-longitude spectrogram 
for Voyager 1 and 2 data (sub-Io system-III longitude of 260°) is shown in 
Figure 1 . This figure displays voyager 1 and 2 data transformed into a 
coordinate system where * >a sub-Io system-III longitude was equal to 260° 

(±3°) The data displayed in Figure 1 is only a subset of the entire Voyager 
data, it is important to note that a band of emission extending in frequency 
from about 6 to 1 3 MHz seen at all Jovian system-III longitudes in Figure 1 
(as well as in Figures 2-4) is partially due to an increase in the sensitivity 
of the PRA instrument at these frequencies (Schauble and Carr, 1983). 

Clearly seen in Figure 1 are two rather narrow emission zones. One 
located at Jovian system-III longitudes between approximately 155°-185° and 


another located between about 305°-330° 


Both of these emission zones extend 


to over 30 MHz. In the region between these two narrow emission zones which 
will hereafter be referred to as northern hemisphere emissions (NHE) is a 
region of les3 distinct but nevertheless structured emission which extends 
between about 2l5°-305° and in frequency from about 5-25 MHz. There is a 
narrow emission feature centered at about 60° and extending in frequency to 
possibly over 30 MHi. but careful examination of this feature reveals that it 
is due to spacecraft noise (interference from other instruments). In 
addition, the emissions seen centered around 140° and extending to about 26 
MHz, and the emissions at about 274° and extending in frequency to almost 40 
MHz, are also of spacecraft origin. 

Also shown on Figures 1-4 are the results of ray tracing emission from 
both northern (circles) and southern (triangles and squares) hemisphere 
sources. The ray tracing results for the southern hemisphere sources are 
shown for initial wave normal angles of 95° (triangles) and 85° (squares) in 
order to investigate the dependence of the resulting arcs on this parameter. 
The agreement for the northern hemisphere emissions is relatively good with 
respect to Jovian longitude. Seen in Figure 1 are model emissions within 
Jovian longitude range 180°-186° and emissions at about 335°. Both of the 
model northern emission cone edges are quite narrow in longitude extent. 

These northern hemisphere model emission cones lie close to the range of the 
observed northern emission cones shown in the same figure, but shifted 
possibility 10° to higher longitudes. The frequency extent for the model 
emission cones is less than observed, however, by about 15 MHz. The reason 
for this discrepancy may to be due to inadequacies in the magnetic field 
model. An increase in the magnetic field strength at the source point from 
about 7 gauss to over 10 gauss would be required for emission at 30 MHz to be 



allowed. This large increase in the field strength suggests either errors in 
the model at low altitudes or transient affects. Since Figure 1 represents 
data taken over a long period of time by two satellites it is possible that 
the observed emissions at the highest frequencies are temporal effects of the 
nonqulescent Jovian magnetic field. 

The model southern emission cones are also shown at specific frequencies 
on Figure 1 superimposed as triangles and squares for clarity. While the 
agreement is not perfect, it is clear that the emissions between the northern 
emission cones shown in Figure 1 are southern hemisphere emissions. The 
frequency extent, is less than that of the NHE and the characteristic 
"arch-shape" is easily discernable. The "left" edge of the arch formed by the 
observed emission from the southern hemisphere is bracketed lay the model 
emissions for wave normal angles of 95° and 85°. The "right" edge of the 
model emission cone is displaced from the observed right edge of the "arch" by 
possibly 30° or so. Changing the initial wave normal angle could produce a 
match with the observations, but so would the introduction of a magnetic field 
model with a different azimuthal dependence. 

Figure 2 is the same format as Figure 1 but now for a fixed sub-Io 
longitude of 160°, in this figure the two observed NHE zones are between 
about 40 0 -90° longitude and between about 195°-240°. The NHE zones are not as 
distinct as they are in Figure 1. There are, for example, reasonably intense 
emissions in the longitude range of 40°-70° and in the frequency range of 
about 17-20 MHz. Another emission reqion appears to be located in the 
longitude range of about 195°-230° and frequency range 16-23 MHz. Theca 
emissions broaden the NHE zones and may represent statistical variations of 



the NHE zones or they may represent emissions not modelled by the present 
assumptions. The results of the ray tracing are shown on Figure 2 as circles 
for the model NHE zones. The model fit is good, but it is apparent that 
the observed "right" NHE zone {at about 235°) is smaller in frequency extant 
than the observed "left" NHE zone (at about 80°), and smaller in frequency 
extent than the model NHE zone. This phenomena, one side of the emission cone 
having a smaller frequency extent than the other, is observed frequently in 
the data and is thus far not adequately modelled. One possible explanation is 
the refractive effects of the Jovian ionosphere which, because of geometry, 
are more extensive on one of the emission cones . This effect is currently 
being investigated. The emissions between the observed NHE cones are not as 
structured as in the example of Figure but it is clear that intense 
radiation persists close to the emission cones produced by the model ray 
tracing and shown as the triangles and circles on Figure 2. While not 
perfect, the qualitative agreement leaves little doubt that the emissions 
between the NHE zones are from the southern hemisphere. For this example the 
right edge of the emission cone is more distinct and appears to be bracketed 
by the model ray tracing results for the two initial wave normal angles. 

The sorted Voyager data shown in Figure 3 is for sub~Io longitude of 
170°. For this example the southern hemisphere emission (SHE) occurs between 
about 95°-21 0° at frequencies less than 25 MHz and has the characteristic 
"arch-shape" described above. The ray tracing results match the SHE very 
nicely, with the "right" edge of the emission cone agreeing better with the 
model for ¥ 0 =85° . The NHE zones are broad and not distinct. The "right" NHE 
zone extends from about 190° to about 310° with a peak frequency occurring 
between 245°-265°. The model ray tracing results fall close to the center of 
this range around 260°. The "left" edge of the observed NHE zone is also not 



very distinct, but appears to be between 35°~80° longitude, with a frequency 
peak possibly between 50°-80°. The model ray tracing results for the "left" 
edge of the NHE zone appear to be shifted by about 15® to 20® to the right 
(higher longitudes) from the observed NHE zone. The data for this sub-Io 
longitude contains a great deal of s/c noise which confuses the match with the 
NHE zones. We have chosen Figure 3 with the confusion of emissions and noise 
to point out the distinctive southern hemisphere emissions. 

The sorted data for sub-Io longitude 150® is shown in Figure 4. The most 
distinctive feature is the band of emission extending in frequency from about 
20 MHz to about 5 MHz and extending in longitude from about 160® to about 
220®. This is identified as the "right" edge of the SHE emission and the 
model ray tracing results bracket this feature reasonably well. The "left" 
edge of the SHE is not as distinct and may overlap with the NHE sons at about, 
75®. The observed NHE zones are not clearly defined. The "left" zone may 
extend from about 35° to almost 100®, and in frequency to over 35 MHz. The 
other NHE zone is even less well defined but may extend from about 200® to 
over 300° with a strong feature at 252®, but emission over 25 MHz is not well 
defined. The model ray tracing results for the NHE at about 75® and at about 
235® are within the observed regions of NHE, but the ray tracing results 
extend to 35 MHz. We have included this figure also to demonstrate the 
distinctive SHE emission even though the NHE zones are much less discrete. 

DISCUSSION 

To our knowledge the results presented here represent the first 
identification of Jovian southern hemisphere DAM emissions. This 
identification has been made possible not only by the comprehensive 



three-dimensional ray tracing utilizing realistic magnetospheric and 
ionospheric plasma and magnetic field models, but also the sorting of the 
Voyager data set with respect to constant sub-Io longitudes. This has allowed 
the quick analysis of many months of Voyager data in a reference frame best 
adapted for study of Io-dependent DAM emissions. The southern hemisphere DAM 
emissions observed are in general lower in frequency and the emission cones 
are narrower. The first of these observations can be explained by noting that 
the southern hemisphere emission sources occur on weaker magnetic field lines 
than northern sources (cf. Acuna and Ness, 1976), and the emission is 
restricted to be slightly Doppler-shifted above the RX cutoff frequency. A 
narrower emission cone results because the Voyager spacecraft, being confined 
approximately to the ecliptic plane, intersects die ■'•■mission cones from the 
northern hemisphere sources closer to the central radius. In other words, the 
emission cones of the southern hemisphere sources are often inside the NHE 
cones. Since the magnetic field is weaker in the southern hemisphere, the 
ratio, fp/fg, will be larger for the same source frequency resulting in 
additional refraction and a smaller emission cone angle. The latter fact also 
explains the "arch-shape" or tendency for the highest frequency emission to 
occur at the same Jovian longitude which appears frequently among the southern 
hemisphere sources. For consecutive frequencies the SHE cones will be 
confined to a smaller range of longitudes than the corresponding NHE cones at 
consecutive frequencies, and the net effect is to produce a bunching of the 
higher frequency ray paths . 



SUMMARY 


we have presented a comparison of the results of three-dimensional ray 
tracing of Jovian DAM emissions from both southern and northern hemisphere 
source positions. The ray tracing is based on a minimum number of assumptions 
which include RX mode emission from source points slightly Doppler-shifted 
above the RX cutoff. A constant initial wave normal angle of 85° was assumed 
for the northern hemisphere sources and both 85® and 95° for the southern 
hemisphere sources, consistent with the results of Staelin, 1981 for electrons 
with energies of several keV. The plasma model is an empirical fit which 
includes the effects of the ionosphere, lo torus, and magnetosphere. The 
magnetic field is the 0-4 model (Acuna and Ness, 1976). The Voyager 1 and 2 
data have been sorted to hold the sub-io longitude constant in order to 
clearly display the lo-dependent features. When directly compared to the 
sorted Voyager data, the ray tracing results clearly distinguish the northern 
and southern emission features, it is found that the southern lo-dependent 
DAM emissions are, in general, lower in frequency, and the emission cone edges 
are observed closer together in longitude along the spacecraft trajectory as 
expected from the model. 

This study shows that the identification of southern hemisphere DAM 
emission establishes strong support that Jovian DAM emission features observed 
by the PRA experiment are primarily due to propagation effects. A simple 
model of wave propagation at large wave normal angles from sources along the 
Io flux tube agrees quite well with many of the Voyager observations. There 
are, of course, discrepancies that need to be addressed such as the observed 
higher upper frequency cutoff of the NHE on Figure 1 , for example . in 
addition, there are other interesting lo-dependent features observed such as 


emission observed at the sub-io longitude of observation (the current lo flux 
tube). We are aware of the power of sorting the data in displaying what would 
otherwise be subtle features in the data set. Attempts to model specific arc 
features can be very frustrating on an individual basis, and the sorting 
technique allows rapid display of qualitative features. We are currently 
introducing new plasma and magnetic field models in an effort to improve the 
model and we will be examining other frequency rangeg by the same procedure. 
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FIGURE CAPTIONS 


Figure 1 . voyager 1 PRA high-band (1,2 to 40.6 MHz) data (from 
February 21 to March 18, 1979) and Voyager 2 PRA high-band data from (June 26 
to July 23, 1979) have been placed in bins of spacecraft system III longitude 
(1965) when lo is within ±3° around 260° system III longitude. Superimposed 
on the voyager data are the model ray tracing results for source points in the 
northern hemisphere (blue circles; ^"SS®) and for source points in the 
southern hemisphere (green triangles ('Jf 0 = , 95°] and red squares [¥ 0 “85®)). Note 
the narrow "arch" of emission between longitudes of about 200° and 300° which 
we believe is due to southern hemisphere sources. 

Figure 2 . The format is the same as Figure 1 but now lo is within ±3® 
around 160° system III longitude. In this figure the "arch" of emission is 
not as distinct as in Figure 1 but there is a "ramp" of emission corresponding 
to one edge of the southern hemisphere emission cone which is bracketed by the 
model ray tracing results between longitudes about 150° and 200°. 

Figure 3 . The format is the same as Figure 1 but with lo centered around 
170° system III longitude. While much of s/c instrument interference is 
present in this figure, the southern hemisphere emission is reasonably 
distinct in its characteristic "arch" between longitudes -of about 100® and 
200° and bracketed by the model ray tracing results . 

Figure 4 . The format is the same as Figure 1 but with lo centered around 
150° system III longitude. Only one edge of the southern hemisphere emission 
cone is distinct between longitudes about 150° to about 200® and bracketed by 
the model ray tracing results . 
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Abstract. A three-dimensional ray tracing coda which incorporates the 0-4 


magnetic field modal (Acuna and Ness, 1979) and a realistic plasma model has 
been used to modal high-curvature decamotric arcs observed by the Voyager 
Planetary Radio Astronomy (PRA) instrument. Two examples of intense, 
isolated, vertex-late, high-curvature arcs were singled out for study. The 
source point wave normal angle was incremented at each of a full range of 
frequencies until the modal rays identically matched the observed arcs of the 
PRA data, several different Doppler shifts were assumed at the source point. 
By this procedure, an accurate relationship between the wave normal angle, 
and the frequency was obtained, with the variation being 70° < ¥ < 80° , and 
^MAX =80° for the arcs considered. As the assumed Doppler shift at the source 
point increased, the value o'i was ^ ounc ^ to decrease, but the general 

shape of the ¥(f) curve was unaffected. By comparing the ray tracing results 
with the independent results of Staelin (1981) regarding emission from beaming 
electrons, it was found that larger Doppler shifts (f/fg > 1*1) at the source 
point produce > 70°, in agreement with the ray tracing results for the 

two arcs considered, only for V|j > 0,7 and Vy/c > 0.32 (where Vy = electron 
velocity parallel to B) . Actual values of Vy/v are unknown, but independent 
observations indicate that vy/c « 0.1. Since the low-curvature arcs are 
believed to result from larger wave normal angles, our results indicate an 
upper limit to the Doppler shift, of Jovian DAM emissions of f/fg < 1.1. 


Introduction 


Frequency-time spectrograms of the voyager I and II low-frequency (1 MHz to 
40 MHz) observations clearly reveal Decametric (DAM) "arc" signatures with 
various radii of curvature (Boischot, 1981} Lecacheux, 1981), In the top 
panel of Figure 1 are displayed DAM arcs observed by Voyager 1 on day 197 of 
1979. The low-curvature arcs seen from about 0510 0T to 0550 UT represent one 
class of event while high-curvature arcs seem to represent a different class. 

A typical example of a high-curvature arc observed by Voyager 1 on day 72 of 
1979 is displayed in the lower panel of Figure 1, extending at both high and 
low frequencies from about 1150 UT to the arc vertex at about 1235 UT. 

Several authors have performed ray tracing analyses utilizing model Jovian 
magnetospheres in order to reproduce the DAM arc features (Goldstein and 
Thieman (1981), Menietti et al. (1984), Green (1984)). These initial efforts 
displayed some success in reproducing the characteristic arc shape with a 
minimum number of assumptions about the emission mechanism. Goldstein and 
Thieman (1981) have found, using an elementary ray tracing analysis, that 
varying the wave normal angle was essential to accurately compare model- 
generated to observed DAM arcs. Goldstein and Thieman obtained an empirical 
relationship to describe the dependence of wave normal angle, on emission 
frequency, f, Their results indicated that ¥ is not a linear function of f 
but rather has a variable peak at a mid-range frequency. These results were 
based on a two-dimensional ray tracing code and represent an initial analysis. 
Hashimoto and Goldstein (1983) published the firBt three-dimensional Jovian 
ray tracing results which offered a convincing explanation of the DAM 
assymmetry of the source position relative to the central meridian longitude 
of lo. A major contribution of their work was in pointing out the 
significance of the ionosphere in refracting DAM emissions. Their analysis, 





2 


which assumed R-X mode emission with sources located just above the R-X cutoff 
frequency (f RX ), indicated that the wave normal angle at the source point 
decreased with increasing ratio of f/f RX at the source point, Hashimoto and 
Goldstein match the occurrence probability of DAM emission in the lo-CML phase 
plana with three-dimensional ray tracing results for three different initial 
wave normal angles. Their results indicate that wave normal angle, 'H - 85° at 
the source point produced the best agreement with ground-based observati ?vf . . 
Menietti et al., 1984 have performed three-dimensional ray tracing of DAM 
emissions and directly compared their results to Voyager observations . They 
point out that for a source region located just above the R-X cutoff frequency 
and for an initial wave normal angle close to 90°, there is good qualitative 
agreement with Voyager data. It appears (e.g. Goldstein and Thieman (1981) and 
Menietti et al,, 1984) that high-curvature arcs such as those shown in Figure 1 
cannot be reproduced with the initial wave normal angle equal to the same 
constant value at all frequencies, even when the refractive effects of the 
Jovian ionosphere are included. The ionospheric model introduced by Hashimoto 
and Goldstein (1983) together with the Acuna and Ness 0-4 magnetic field model 
will affect northern hemisphere sources at mid to high latitudes only for f > 

20 MHz. ionospheric refraction cannot explain the low frequency refraction 
necessary to produce the observed signature of the high-curvature arcs. 

In this paper, three-dimensional ray tracing is used to examine variable 
wave normal angle effects on the resulting DAM arc structure at frequencies 
just above the R-X cutoff and with several assumed ratios of f/f g at the source 
point. We use whatever wave normal angle is necessary to match Voyager 
observations of two intense, vertex-late, high-curvature DAM arcs. The results 
indicate that one possible explanation of the double-valued frequency nature of 
the high-curvature arcs is Doppler-shifted gyroemission from a beam of 


electrons with E » 10 keV. 


An independent investigation of Doppler-shifted gyroemission utilizing the 
results of Staelin (1981) has allowed us to set an approximate upper limit on 
the Doppler shift of the stimulating electrons. Wave normal angles which fall 
in the range of the ray tracing results, i,e< ^MAX > 70° are possible for 
Doppler-shifted gyroemission for f/fg > 1.1 only if V n/V > .7 and Vy/c > .32. 
This work represents the first comprehensive three-dimensional ray tracing 
investigation of varying the wave normal angle and Doppler shift utilizing 
realistic Jovian plasma and magnetic field models with direct comparison of the 
results to Voyager observations of DAM arcs. It thus complements and extends 
the results of Goldstein and Thieman (1981) and the work of Hashimoto and 
Goldstein (1983). 

The ray tracing code used in this study has been described by Menietti et. 
al,, (1984) so only a brief description will be given here. The code is 
three-dimensional snd based on the cold plasma formulation of stix (1962) and 
integrates the Haselgrove (1955) equations. The magnetic field model is the 
0-4 model (Acuna and Ness, 1976) and the plasma model is that introduced by 
Sentman and Goertz (1977). The three-dimensional results reduce to the 
two-dimensional results of Green et. al. (1977), Green and Gurnett (1980), and 
Lecacheux (1981). 

The assumed emission mechanism requires that: (1) the radiation is R-X mode; 

(2) the source region is at the foot of the instantaneous io flux tube and the 
frequency of the emission in the source region is greater than the Rx cutoff 
frequency, and (3) the emission cone is hollow. 

High-Curvature Arcs 

Often seen in the Voyager PRA data are high-curvature arcs extending over 
time intervals of thirty minutes or more with frequency ranges of less than 1 
MHz to over 30 MHz, (Figure 1 for example). These arcs occur isolated and 


nested in groups. Two examples of isolated high-curvature arcs are shown 
plotted in Figure 2. These were observed by Voyager 1 on day 72 of 1979 and by 
Voyager 2 on day 188 of 1979. In order to model this class of arc the wave 
normal angle, f, (the angle between 1c and the magnetic field, B) was varied at 
the source point for each chosen frequency. Source points were chosen at a set 
of frequencies: 2, 5, 10, 15, 20, 25, and 30 MHz from the condition f/fq = 1.02 
and f/fg 3 1.2. These ratios were chosen to be consistent with Doppler-shifted 
emission at large wave normal angles. Gyroemission is not assumed for the ray 
tracing. The condition f > f KX was satisfied at each source point. For a 
specific frequency, ratio f/f g , and sub-Io longitude (lo flux tube) an initial 
angle, ¥, was chosen to begin ray tracing and the proximity of the resulting ray 
path to the spacecraft position (already known from observation of a particular 
high-curvature arc) was noted. If the resulting ray path did not intersect the 
spacecraft position, a new initial wave normal angle, ¥, was chosen and the 
process repeated. A new io flux tube was chosen for each frequency in 
agreement with observations because lo revolves slightly relative to Jupiter 
during the time of observation of an arc. The longitude of the currently 
stimulated lo flux tube was chosen according to observations. Since the 
observed arcs were treated as input constraints to the ray tracing process, 
this procedure generated model arcs which exactly matched the observed arcs 
when ¥ varied with frequency according to the curves shown in Figures 3 and 4. 

Results 

, The results for the two high-curvature arcs observed by Voyager 1 on day 72 
of 1979 and on day 188 of 1979 are shown in Figures 3 and 4 respectively. As 
can be seen, the ¥(f) curves for each modelled arc are qualitatively similar 
(peaking at f » 10 MHz and falling off at lower and higher frequencies). The 
curves of Figure 3 were evaluated for f/f g = 


1.02 and 1.2 at the source 
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points. An increase in the Doppler shift from f/f g = 1.02 to f/f g » 1.2 
shifts the curve downward approximately 3° while the characteristic shape is 
maintained. The data for day 188 in Figure 4 for f/fg “ 1*2 and f/f g = 1«02 
show the same general functional relationship between ¥ and f as shown in 
Figure 3. In Figure 4 the shift in the curves is approximately 8° at the 
highest frequency. The arc vertex frequency (~ 10 MHz) requires the largest 
wave normal angle (~ 79°) at f/f g “ 1.02 as determined for both day 72 and day 
188. 


For both high-curvature arcs plotted in Figure 2 the lo flux tube 
intersects Jupiter at a longitude where the azimuthal variation in the B-field 
(0-4 model) is slight. The results may be explained by noting that since the 
source positions for the highest frequencies along the lo flux tube are closer 
to the planet by about 1 Rj, increased refraction of the higher frequencies 
refracts the emission cones so that the spacecraft observes the higher 
frequencies "inside” the emission cones of the lower frequencies (cf. 

Goldstein and Thleman, 1981). 

Doppler-Shif t Gyroemission Discussion 

Let us consider the case of Doppler-shifted gyroemission from high energy, 
field-aligned electrons. The ray tracing analysis is not dependent on an 
emission mechanism and so this section is distinct from the ray tracing 
results. Staelin (1981) has used the results of Goldreich and Lynden-Bell 
(1969) to obtain an expression for the maximum wave normal angle expected for 
Doppler-shifted R-X emission (assuming the emission phase velocity is 
approximately equal to c) as 


¥ 

MAX 


co,-’ i ligt/alinl ) 


( 1 ) 
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where K = Vy/V and Y is the relativistic term. The Doppler-shifted frequency 
is given by 

= 1 + cos (2) 

fg C 2 C 

For a given ratio, Vy/c, the above equation 2 yields the Doppler shift ratio, 
f/fg. Then, for a given value of K, y is determined and can be found. The 

precise values of K and Vy appropriate for Jupiter are unknown. Using Vy «> .1c 
(as Staelin suggests is consistent with observed de came trie S burst frequency 

drifts) and choosing an intermediate value of K » .5 gives y = 1.02 and “ 

78° (if we assume fp/f g << 1 which is the case near most of the source points 
considered) . This value of 'J'^x c ^ ose to the value of the peak wave normal 
angle shown in Figures 3 and 4 for f/fg ~ 1.02. Figure 5 contains a family of 

^MAX versus K curves for values of Vy/c which vary from 0.1 to 0.45. This 

figure differs from Figure 2 of Staelin (1981) in that we have labelled curves 
as a function of the assumed Doppler shift ratio as shown in Table I. The value 
of y varies slightly for each curve as shown. These curves demonstrate that the 
larger Doppler shifts (as shown in curves 4 and 5 for f/fg = 1*1 and f/fg = 1.2 
respectively) produce values of ^m^x > 70 ° consistent with the ray tracing 
results shown in the curves of Figures 3 and 4 only if K > .7 and Vy > .32c. 
Accepting Vy a .1c as a reasonable value, our results indicate that f/f g ** 1.1 
is an upper limit to the Doppler shift at the source point for the arcs 
considered in this study. 


TABLE 1 


PARAMETERS USED IN FIGURE 5 

Curve 

f/fg 

V,|/c 

Ymin 

YMAX 

1 

1 .01 

0.10 

1 .01 

1 .03 

2 

1 .02 

0.14 

1 .01 

1 .04 

3 

1 .04 

0.20 

1 .02 

1 .09 

4 

1 .10 

0.32 

1 .05 

1 .12 

5 

1 .20 

0.45 

1.10 

1 .20 



Conclusions 




In the Voyager PRA data set one frequently observes either isolated or nested 
groups of arcs that have a larger curvature and are difficult to explain if one 
assumes Hf => constant « 90°, Two examples of isolated, vertex-late, 
high-curvature arcs have been matched through ray tracing by allowing the wave 
normal angle to vary with frequency. The Uf(f) relationship that produced these 
spectral features was peaked around 10 MHz and fell off on either side in 
qualitative agreement with the earlier two-dimensional ray tracing results of 
Goldstein and Thieman (1901). 

If a gyroemission mechanism is assumed that produces Doppler-shifted 
emission, the results of Staelin (1981) predict a frequency dependence of the 
wave normal angle. An energetic, field-aligned electron beam (V|| = 0.1c) could 
produce maximum wave normal angles in agreement with those found in the present 
models if f/f g « 1.02 and K « 0.5. Such energies (10 keV) are comparable to 
suggested energies of several keV for accelerated electrons near the surface of 
lo which may be the source of DAM emission (cf. shawhan, 1976; Goldstein and 
Goertz, 1983). We emphasize, however, that our ray tracing results are not 
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dependant on a gyroemission mechanism. Each of the 'F(f) curves were obtained 
from ray tracing by requiring that the model DAM high-curvature arcs matched 
those observed by Voyager 1 on day 72 (Figure 3) and day 188 (Figure 4) of 1979. 
The characteristic shape of the curves (with a peak at the vertex frequency) 
does not significantly change as the Doppler shift increases. The wave normal 
angle, Y, varies between 70° and 80° for all Doppler shifts considered with 
smaller values of ¥ resulting from larger assumed Doppler shifts. Independent 
of the ray tracing results, the condition that gyroemission occur for f > f RX 
(equation 1 [cf. staelin, 1981}) suggests that an electron beam with K » .5 and 
V || “ .1c will result in which agrees favorably with the results shown in 

Figures 3 and 4. Figure 5 and the values listed in Table 1 indicate that larger 
Doppler shifts (f/fg 3 1.1 and f/fg a 1.2) will produce > 70° in agreement 

with the ray tracing results only for K > .7 and Vy/c > .35. The valuo of K = 
Vy/V appropriate for Jupiter is not known. But observations of frequency drifts 
of Jovian S bursts imply that Vy 3 .1c. As seen from Table 1 this implies that 
for the two hiqh-curvature arcs presently modelled, curves 1 and 2 for Vy/c < 

,14 and f/fg 3 1.02 best match the ray tracing results. The two arcs presently 
modelled were chosen because of the large frequency dispersion and intense 
emission, we believe these arcs represent emission at close to the smallest 
wave normal angles and the largest Doppler shifts for observed DAM arcs 
(Figures 3 and 4 indicate that a larger assumed Doppler shift at the source 
point yields a smaller value of • 0ur results thus indicate that f/f g « 

1,1 is the largest to be expected for Jovian DAM arcs, with typical values being 
f/fg < 1.02. While we have chosen only vertex-late arcs for this investigation 
(due to the intensity of the emission and large frequency dispersion) we expect 
no differences in our conclusions for the vertex-early arcs . We have excluded 
radiation from Alfven wave reflections located away from the Io flux tube. If 
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simultaneous emission from multiple flux tubes was responsible for the 
production of the two isolated arcs modelled, our conclusions would have to be 
modified. 
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Figure Captions 


Fig. 1. Frequency-time spectrograms of Voyager low-frequency radio 
emissions. The emission intensity is gray-shade coded acording to the bar 
graph at the left. In the top panel low-curvature DAM arcs observed by Voyager 
1 on day 197 of 1979 are seen from about 0510 UT to about 0550 UT. in the 

bottom panel is seen a typical example of a high-curvature DAM arc obsrerved by 

Voyager 1 on day 72 of 1979 and extending from about 1150 UT to arc vertex at 
about 1235 UT. 

Fig. 2. A frequency versus system-III longitude plot of Voyager data for 
two observed isolated DAM high-curvature arcs: Voyager 1 for day 72 of 1979 
(shown in spectrogram in Figure 1); and voyager 2 for day 188 of 1979. 

Fig. 3. Initial wave normal angle, versus frequency obtained from ray 
tracing which resulted in a DAM arc identical in curvature to that observed by 
Voyager 1 on day 72 of 1979. The upper curve is for an initial Doppler-shift 
of f/fg = 1.02, while the lower curve is for f/fg = t»2« 

Fig. 4. same format as Figure 3 except for comparison to the arc observed 

by Voyager 2 on day 188 of 1979. 

Fig. 5. Maximum wave normal angle, versus K ( V |j/V) for a family of 

curves with parameters varying as shown in Table 1 . The ratio, f/fg , increases 
from 1.01 (curve 1) to 1.2 (curve 5). Gamma is approximately a constant for 


each curve. 
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ABSTRACT 


The post-encounter Voyager detection of a night-side deca- 
metric emission which differs from pre-encounter and earth-based 
observations is examined in this study. This intense radiation 
is observed in the frequency range from about 12 to 25 MHz, is 
detected in bursts lasting a couple of hours, is localized in 
Jovian longitude, and is clearly an Io phase independent emis- 
sion. A spatial survey of this emission reveals that it is 
observed when the Voyager spacecraft are in the System III 
(1965) longitude range from 150 to nearly 190 degrees, at posi- 
tive magnetic latitudes, and with a local time from 0 to 3.9 
hours on the night-side of the planet. The night-side nature of 
this emission prevented its detection by Earth-based radio 
observations which are made in the local time range from about 
11.3 to 12.7 hours. 

This night-side component of decametric radiation may be 
coming from the low-altitude night-side auroral zone which is 
linked to the Jovian plasmasheet by magnetic field lines that 
are at higer latitudes than Io' s field lines. From all indi- 
cations, it appears that this emission is the Jovian counterpart 
to the Earth's auroral kilometric radiation or AKR, hence it is 
referred to here are Jovian Auroral Decametric Radiation or 
JADR. Since JADR is seen nearly every time the Voyager space- 
craft are in the emission pattern of the radiation, the night- 
side Jovian tail field auroral zone may be continuously active. 
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INTRODUCTION 


It is cleat from satellite and Earth-based observations 
that the Jovian magnetosphere provides a rich environment for 
the generation of many types of free escaping radio emissions. 
Earth-based radio measurements as early as 1955 discovered 
emission from Jupiter at decametric or DAM wavelengths. 

Brigg discovered in 1964 that the phase of Io has a strong 
influence on whether or not decametric radio emissions from 
Jupiter are detected at earth. This Io effect has led to the 
classification of the decametric emissions as Io dependent and 
Io independent. The Io independent emissions are less intense 
and show Jovian longitude control. 

This paper will discuss observations from the PRA experi- 
ment onboard the Voyager 1 and 2 spacecraft of a new type of Io 
independent emission at decametric wavelengths which shows a 
strong dependence on the observing longitude and magnetic lati- 
tude and are seen only on the night-side of Jupiter. 
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OBSERVATIONS 


In order to concentrate on the Io independent emissions 
close examination of the data must be made to clearly identify 
the Io dependent emissions. Decametric measurements have been 
routinely put into a coordinate system of central meridian 
longitude (CML) versus Io phase (see for example Alexander et 
al., 1981). In this coordinate system Io dependent emissions 
have been classified as "sources" or regions of CML and Io phase 
where the probability of measuring decametric radiation is high. 

It has been known for sometime that the most intense decametric 
emissions are from sources A, B, and C. Green, 1984 using 
Voyager 1 and 2 observations over several months organized the 
Voyager data into spacecraf t-Io longitude difference versus ■ 

frequency for fixed positions of io in System III (’65) longi- 

tute. As discussed in detail in Green, 1984 the purpose of 
plotting the Voyager data in this new coordinate system is to 

J 

eliminate the relative motion of Jupiter and Io. The result of j 

. ' ! 

the coordinate transformation of the Voyager data (see Figures | 

3, 4, 5, and 6 of Green, 1984) is the same as would be observed jj 

by a spacecraft that . moved completely around Jupiter at a con- 
stant radial distance keeping io and Jupiter fixed. It is 
important to note that out of several months of observations 
only a relatively small amount of data is used which satisfied j 
the conditions of having Io at a specified system III longitude. 
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The analysis presented here carries the work of Green, 1984 
a step further by selecting a band of frequencies for each fixed 
Io sorted data plot of spacecraf t-Io longitude difference (or 
phase) (in this case from 15 to 21 MHz). Figure 1 presents the 
occurrence of decametric emission in the frequency range 15 to 
21 MHz over the period June 26 to July 23, 1979 versus Io longi- 
tude. The top panel is determined from all PRA observations on 

2 

the inbound and outbound leg of the Voyager 2 flyby. A 1/R 

factor was used to take into account the radial dependence of 

the emission. An occurrence block is shown when the average 

decametric intensity (normalized to 10 ) over the 15 to 21 MHz 

J 

-17 2 

frequency range was greater than 1.1 x 10 Watts/ m /Hz. 
Observations closer than 30 Rj to the center of Jupiter were not 
used to avoid any possible propagation effects in the inner 
y* Jovian magnetosphere. All spacecraft observations are organized 

in longitude bins 3 degrees wide and plotted such that lo’s 
position is along a vertical line at the center of the plot, 
with the difference between spacecraft longitude and Io longi- 
tude increasing to the right of center and decreasing to the 

left of center. Along the ordinate the data is organized into 6 
degree wide Io longitude increments spaced every 10 degrees 
apart. Labeled in all three panels of Figure 1 are the expected 
locations of the Io dependent sources A, B, and C. Vertical 
features on these plots would maintain constant Io-observer 
phase and hence, be labeled "Io phase dependent". Note that 
periodic observations of decametric radiation are seen extending 


from the upper left corner of Figure 1 (top panel) to the lower : 
right. Although these observations cross into the B source it 
is certain that they are not Io phase dependent. 

The middle and bottom panels of Figure 1 are plotted in the 
same format as the top panel, but contain only inbound data 
(middle panel) and outbound data (bottom panel) data. When the 
Voyager 2 data are separated in this manner, incomplete sampling 
occurs at many of the longitudes. The regions of no observa- 
tions in this coordinate system are shown shaded. The separa- ' 
tion of the Voyager 2 data into inbound and outbound passes is j 
performed in order to examine the local time effects in the 

observation of the periodic decametric radiation pointed out in ; 

\ 

the top panel. The PRA measurements on the inbound portion of jj 

jj 

the Voyager flybys restrict the majority of the observations to j 
the day-side Jovian magnetosphere from 9.5 to 11.0 hours local 
time (LT) and the outbound from 0 to 2.8 hours LT. A comparison j 
between the middle and bottom panels indicates that the periodic | 
emission seen in the top panel extending from upper left to i ^ 

i 

lower right is observed only on the outbound leg of the Voyager 
2 flyby. This analysis confirms the detection of a night-side 
component of decametric emission from Jupiter, which is clearly! 

Io phase independent. 

The spectrogram of Figure 2 shows the frequency- time char- 
acteristics of a typical night-side decametric burst lasting z 
couple of hours. The enhancement in decametric emission at the 
low frequencies, from about 8 to 16 MHz, may be attributed to at 
increase in sensitivity of the PRA instrument in this frequency 
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range recently determined by Schauble and Carr, 1983. The 
enhanced sensitivity in this frequency range sometimes gives the 
appearance that a broad Jovian emission band exits. In order to 
avoid the sensitivity effect and to provide a concise picture, 
the entire frequency range of the night-side decametric emission 
was not used (only frequencies from 15 to 21 MHz were used) in 
the determination of DAM occurrence of Figure 1. Usually in the 
8 to 16 MHz frequency range, "arc" structures (see Warwick et 
al., 1979a and Leblanc, 1981) are found that are typical of 
non-Io emissions, however, in Figure 2, at higher frequencies 
(>12 MHz) this emission is not made up of a series of arcs and 
is much broader (as observed in time) than a single arc. The 
frequency- time structure of this emission is seen periodically 
in the Voyager 2 observations from about 25 R ^ until the emis- 
sion is just barely measured above the receiver's noise level at 
distances greater than 120 R j . 

The spectral characteristics of the night-side decametric 
emission of Figure 2 is shown in Figure 3 where the power flux 
emission (expressed in Watts/m /Hz) has been normalized to 10 
R j . In addition, Figure 3 shows only the maximum power flux 
observed at ea;;h frequency over the spacecraft longitude range 
from 150 to 160 degrees (about 17 minutes of data) . Note that 
this emission reaches peak intensity from 14 to 25 MHz and 
decreases rapidly at higher and lower frequencies. This spec- 
tral shape is not characteristic of the Io dependent decametric 
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emission (Carr et al., 1983) but it does more closely resemble 
the spectral shape of the Earth's auroral kilometric radiation 
(see for example Gurnett, 1974). 

The repetition of the night-side decametric bursts may be 
understood by examining if there is any dependence on the space- 
craft position in magnetic latitude and longitude for observing 
the emission. On a "wiggle" plot (position in magnetic latitude 
and radial distance with time) , of the Voyager 2 outbound tra- 

. i 

jectory in Figure 3, are dark strips where complete night-side! 

j 

decametric bursts are observed. These bursts are observed exclu-j. 
sively at high positive magnetic latitudes with the high fre-l 

>j 

quencies seen at magnetic latitudes greater than about 9 de-jj 

H 

grees. It is important to note that the highest frequencies of-J 

jl 

the night-side decametric emission (20 to 22 MHz) are not mea-j; 
sured when the spacecraft is at the highest magnetic dipolef 
latitude (corresponding to approximately 200 degrees longitude)! 
but at the latitude which corresponds to the longitude (abou ? j 

I 

150 degrees) of the believed maximum surface magnetic fielt:*-„ 
strength (see for example, Connerney, et al . , 1981). This rela 
tionship is very suggestive of low altitude source regions whic 
have a strong dependence on Jupiter's magnetic field. In addi . 
tion, since the maximum frequency is observed at or very nes 
the longitude of the maximum field strength the beaming of this 
emission appears to nearly parallel to the magnetic field in ti 
source region. This situation is very different from the bro. 
emission cones (large initial wave normal angles) that a: 
believed to be so characteristic of the Io A and B sources (s 
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Green/ 1984) . Better understanding of the emission and propaga- 
tion characteristics of this radiation is currently under study 
using ray tracing techniques. 

The spatial distribution of the night-side DAM as suggested 
by Figure 1 can be further investigated by using Voyager 1 PRA 
data as well as the Voyager 2. The Voyager outbound trajec- 
tories, looking down from the north onto an equatorial plane, 
are plotted in Figure 5. The shaded regions are at the position 
of the spacecraft where the night-side DAM was measured. The 
cross hatched block at the expected location of the emission 

along the Voyager 2 trajectory is due to the inability to clear- ] 

j 

ly separate the night-side DAM emission from intense Io depen- \ 

j; 

dent source B radiation measured at the same time. It is impor- I 

I 

. ; j tant to note however, that the Voyager 1 measurements of the j, 

? night-side DAM abruptly stopped after a radial distance of about j 

58 Rj and local time of 3.9 hours. The Voyager 2 observations ] 
of this emission, shown in Figure 5, continue until the inten- ' 
sity of the emission is near the PRA experiment noise level at „ 
radial distances greater than 120 R.. The overall repetition of 

3 j 

the Voyager 2 observations (and to some extent the Voyager 1) j 
strongly indicates that the absence of Voyager 1 night-side DAM) 
observations past 58 Rj on the outbound trajectory is due to a) 
spatial and not temporal effect and that the this emission has s| 

distinct outer boundary that Voyager 1 passed through. Combined 

i ■ 

with the results of Figure 4, the "inner boundary" labeled ir 

i 

Figure 5 is probably due to the low magnetic latitude measure-] 
ments near the encounter. 1 


DISCUSSION 


The night-side DAM emission reported here has been com- 
mented on in the paper by Alexander et al., 1981. This emis- 
sion, as described here, can be easily seen in the synoptic 
analysis of Plate 1 in the Alexander et al., 1981 paper for only 
the Voyager 2 after encounter plot with traces of the emission 
seen in the Voyager 1 after encounter data. The low occurrence 
probably of the night-side DAM from the Voyager 1 observations 
in the Alexander et al., 1981 paper is undoubtably due to the j 

f 

fact that Voyager 1 was in the emission region for only a short j; 

i) 

time (see Figure 4 of this paper) . Alexander et al., 1981 did •; 

ij 

point out that the " Io- independent emission at frequencies above | 
20 MHz exhibit marked local time effects". Their interpretation 
of this emission as a feature which "probably corresponds to the 
weak Io- independent B source recorded in earth-based surveys" is 
not the interpretation presented here. It is without question 
that the night-side DAM can only be measured at high magnetic f 
latitudes (greater than 9 degrees) and on the night-side ofj 
Jupiter, both of which, the Earth-based measurements are unable 1 
to achieve. The situation with the lower frequency extension of 
this emission (less than 12 MHz) however is not at all clear and 
may very well be related to the Io- independent B source as 
discussed by Alexander et al., 1981. Additional work needs tc 
be done to clarify this point. 


From the above observations it is important to compare the 
night-side DAM to other well known magnetospher ic emissions in 
order to gain insight on its possible origin and generation. It 
is generally believed that the earth's auroral kilometric radia- 
tion (AKR) is the terrestrial counterpart to the decametric 
radiation from Jupiter. These emissions have many physical 
similarities such as the same polarization, the same frequency 
range when normalized with respect to the planets magnetic 
field, high latitude low altitude source regions, and the most 
intense free escaping radio frequency radiation generated in 
their respective magnetospheres. Many generation mechanisms 
have been proposed for both AKR and DAM which are the same (see 
for example Melrose, 1976). However, when considering the Io 
control of the intense DAM radiation there is no obvious terres- 
trial counterpart. A more analagous comparison would be between 
AKR and the Io independent emissions. The night-side DAM emis- 
sions, discussed here, are the only Io independent emission 
which have spatial and spectral characteristics similar to AKR. 
The angular distribution of AKR has been extensively studied 
(see Green et al., 1977, for example) and it is clear that AKR 
emission exhibits latitudinal propagation effects and is seen 
primarily on the night-side of the earth. If the night-side DAM 
is the Jovian counterpart to AKR, then it may be coming from 
Jupiter's tail field auroral zone (with low altitude sources) 
which is at higher latitudes than the field lines that Io inter- 
sects (Connerney, et al . , 1981). Because of these similarities, 
the night-side decametric emissions will be referred to here as 
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Jovian Auroral Decametric Radiation or JADR. The persistence of 
JADR in the Voyager 2 data may therefore indicate that Jupiter 
has an extremely active and continuous night-side particle 
precipitation in the Jovian tail-field auroral zone. It is 
important to note that recent work by Maeda and Carr, 1984 
suggests that the Io-independent A source, seen on the day-side 
of the Earth, may have an Jovian auroral zone origin. How this 
emission relates to the JADR discussed here is not known and it 
is hoped that future ray tracing studies may be successful in 
determining the relationship between the Io-independent A source 
and JADR. 
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CONCLUSIONS 


We have presented an identification of one of the deca- 
metric components generated in the Jovian magnetosphere. This 
emission, referred to here as JADR, has been measured by the 
Voyager 1 and 2 PRA experiments and has the following major 
characteristics : 

1. JADR is an Io-phase independent emission in the fre- 
quency range from 12 to 25MHz. 

2. JADR (from 16 to 25MHz) cannot be categorized as 
having an "arc" like structure characteristic of much 
of the lo-dependent and -independent emissions. 

3. The JADR power spectrum is not typical of the 
lo-dependent or Io- independent DAM emissions. 

4. There is a strong latitudinal and longitudinal beam- 

ing of JADR since it is observed at positive magnetic 
latitudes with the highest frequencies measured when 
the largest surface magnetic field strength (near 150 
degrees) is near the longitude meridian of the space- 
craft. ' s < 
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5, It is observed only on the outbound portions of the 
Voyager flyby trajectories. 

6. A strong night-side beaming effect is implied when an 
abrupt halt to intense JADR observations from Voyager 
1 occurs at local times greater than 3.9 hours, while 
Voyager 2, at local times less than 2.8 hours ob- 
serves the emission to radial distances greater than 

120 R j . 


A number of conclusions can be made from the observed 
characteristics of the JADR emissions and their similarity to 
other magnetospher ic emissions. These include: 

1. A possible propagation effect may be responsible for 
the magnetic latitude asymmetry seen. 

2. Any southern hemisphere JADR emissions would only be 
observed at the lowest frequency (12 MHz) since neither 
Voyager 1 or Voyager 2 were at southern magnetic lati- 
tudes greater than about 7 degrees. 

3. JADR emission may be the Jovian counterpart to the 
earth's auroral kilometric radiation and emanate from a 
continuously active Jovian tail field auroral zone. 


<»%. -lAJk. 
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Figure 1 


Occurrence of DAM emission from Voyager 2 spacecraft is 
shown by the blacken squares when the average intensity is 
greater than 1.1 x 10 Watts/m /Hz. {normalized to 10 Rj) 
over the frequency range from 15 to 21 MHz. This three panel 
figure is derived from observations taken along the outbound 
(bottom), inbound (middle), and both inbound and outbound (top) 
Voyager 2 trajectory. The DAM occurrences are plotted in degree 
difference between the spacecraft and Io versus Io System III 
(’65) longitude. In the top and bottom panels only, a systema- 
tic DAM non-Io emission is observed from the upper left to the 
lower right corner of the plot. 
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Figure 2 


This Voyager 2 frequency-time spectrogram shows the basic 
characteristics of the JADR emissions. A night-side decametric 
burst usually lasts about an hour and is first seen in Voyager 1 
and 2 data at high frequencies (typically 21 to 22 MHz) and as a 
function of time to lower and lower frequencies until a broad 
emission band appears from 8 to 12 MHz. 





m 
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Figure 3 

The maximum power flux spectrum of the night-side deca- 
metric emission is shown in this figure. Note that the spectral 
peak occurs from 14 to 22 MhZ and rapidly decreases at higher 
and lower frequencies. 


POWER FLUX (W Hz' 


VOYAGER 2 79/1 94 0825 TO 08A2 SCET 
S/C LONGITUDE 150°-1 60° 



Figure 3 


•> 





23 


Figure 4 


The position of the Voyager 2 spacecraft when the night- 
side DAM emissions are observed in magnetic latitude and radial 
distance coordinates are shown here as darkened strips during 
the outbound portion of its trajectory. The characteristic 
shape of these emissions as pointed out in Figure 2 were used in 
determining the spacecraft positions of the shaded strips. Note 
that the night-side DAM emissions are all seen at positive 
magnetic latitudes with the highest frequencies observed near 
the beginning of each strip. 
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Figure 5 


The dark block regions along the outbound Voyager 1 and 2 
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